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Abstract 

We investigate finite size effects of the pion matrix element of the non- 
singlet, twist-2 operator corresponding to the average momentum of non- 
singlet quark densities. Using the quenched approximation, they come out to 
be surprisingly large when compared to the finite size effects of the pion mass. 
As a consequence, simulations of corresponding nucleon matrix elements could 
be affected by finite size effects even stronger which could lead to serious 
systematic uncertainties in their evaluation. 



1 Introduction 



A complete description of moments of parton distribution functions (PDF) from first 
principles is still missing. In order to provide reliable numerical determinations, it 
is necessary to control all systematic uncertainties which include the continuum 
limit, non-perturbative renormalization, chiral extrapolation, finite size effects and 
quenching. Several studies on hadron matrix elements of twist-2 operators have been 
performed since now 12] • In |2j a study of the continuum limit employing a per- 
turbative renormalization shows small cutoff effects. Applying a non-perturbative 
renormalization [3], slightly reduces an observed discrepancy between theory and ex- 
periment [3]. First computations in the unquenched case jT|, at a rather high value 
of the quark mass, seems to indicate no quantitative differences with the quenched 
case. In (3] it has been suggested that the discrepancy between theory and experi- 
ment could be explained by a rather strong dependence of the matrix elements on 
the quark mass in the chiral region when a very small quark mass is entered. The 
simulations performed since now are carried out at values of the quark mass too 
far away from the physical point to test this anticipated mass dependence. One 
problem that to our knowledge has not been faced in a systematic manner is how 
the finite size of the lattice can influence such computations. 

In this letter we present a study of the finite size effects (FSE) of the quenched 
pion matrix element of the twist-2 operator corresponding to the average momentum 
of non-singlet quark densities. Preliminary results were presented at the last lattice 
conference jH]. 

It has already been noticed some time ago [7] and more recently in [HJ E] , that 
FSE expecially in the unquenched case can be a serious source of systematic error. 
In this letter we want to emphasize that even in the quenched case, FSE for pion 
3-point functions involving twist-2 operators, can be much larger than the FSE for 
the hadron masses. Thus, there is a clear warning that for the unquenched case and 
moreover for nucleon matrix elements such FSE should be studied very carefully in 
order to avoid an underestimation of this systematic effect. 

2 Setup and basic definitions 

The moments of parton distribution functions (PDF) (x^ N ~^} are related to matrix 
elements of leading twist r (r =dim-spin) operators O of given spin, between hadron 
states h(p) 

(h(p)\0^ m \h(p)) = M^in^-.-p^ 

+terms 8^ (1) 

(x^)^)=M^\^ = Q), (2) 



1 



where the renormalization scale /i is identified with the momentum transfer Q in 
deep inelastic scattering (DIS). 

Our setup of lattice QCD is on a hyper-cubic euclidean lattice with spacing a and 
size L 3 x T. We impose periodic boundary conditions in the spatial directions and 
Dirichlet boundary conditions in time, as they are used to formulate the Schrodinger 
functional (SF) (we refer to these references for unexplained notations). In 
this paper we will consider homogeneous boundary conditions, where the spatial 
components of the gauge potentials at the boundaries and also the fermion boundary 
source fields are set to zero. In this case the Schrodinger functional partition function 
can be written as [TTj 

Z = (*o|e-™P|*o> (3) 

where the initial and final states |io) carry the quantum numbers of the vacuum. 
We have studied the matrix elements between charged pion states of the following 
operator 

O u (x) = -$(ar) 74 B 4 ^7fe D k 4>(x) . (4) 
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We indicate with ( and ( (and the corresponding (' and £') a flavour doublet. As 
discussed in jT2| we recall that (, (, (' and (' are functional derivatives with respect 
to the boundary source fields. 

The states with charged pion quantum numbers in the Schrodinger functional 
are the dimensionless fields 

S = J} £ C(y)75T + C(z) S' = ^ E C'(u)75T-C'(v) (5) 

where = -^(r 1 ± ir 2 ) and r k with k = 1, 2, 3 are the usual Pauli matrices. The 
correlation function to extract the matrix element is 

f u (x ) = (SO u (x)S') . (6) 

The Wick contractions of this correlation function contain also a disconnected piece 
that we neglect consistently with the fact that we will do quenched simulations. For 
normalization purposes it is important to define 

fx = -|<ssr> . (7) 

The correlation functions fa and f\ have the following quantum mechanical repre- 
sentations 

fu(x ) = Z- l {i^ T -^ m ¥0^x)e-^ m ¥\i^) (8) 
Zi = ^(v|e- TH P|v>. (9) 
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Inserting a complete set of eigenstates of the Hamiltonian and retaining only the 
first non-leading corrections we have 

/ 4 4(*o)^p 2 (O,7r|0 44 (x)|O,^ (10) 



f x ~ p 2 e- m « T . (11) 

The matrix element we are interested in, neglecting contributions from excited 
states, can then be extracted from the plateau value of the following ratio: 

M^- = (0,n\O u (x)\0,n). (12) 

h 

In order to relate this numerically computed ratio with the corresponding continuum 
operators in the Minkowski space, we need a suitable normalization factor 

(a;) = — (0,7r|(944(z)|0,7r) , (13) 

where we always take for m n the infinite volume pion mass m 7r (L = oo), summa- 
rized in table El in order to disentangle the FSE of the pion mass and (x) itself. 
Furthermore, we consider only the bare matrix elements without renormalization 
factor, since the FSE remain unchanged when employing the appropriate Z factors. 

Using the procedure of ref. JT] we have extracted the plateau values for the 
effective pion mass and the 3-point function with a systematic relative error coming 
from the excited states of 0.1%, respectively 0.4%, well below the statistical accuracy 
of our computations. Technical details can be found in a forthcoming publication 

B3|. 



3 Numerical details and results 

We have performed a set of quenched simulations at fixed lattice spacing a = 0.079 
fm (/3 = 6.1), for three values of the quark mass (ac = 0.1340, 0.1345, 0.1350), and 
several lattice sizes, using the non-perturbatively improved clover action |14j . For 
some simulations we have also used two different time extents of our lattice (T, T', 
with X" < T). Our simulation parameters are summarized in tabled Using r = 0.5 
fm to fix the scale [15] . we have lattice sizes varying between 0.9 fm and 2.5 fm, and 
values of z — m n L varying between 2.5 and 11.2. 

We have made fits of the data both with an exponential and power dependence 
on the lattice size L (in the following equation we generically indicate with F(L) 
both and (x)): 

exp : F(L) =a + — ^ exp{-a 2 L} (14) 
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Figure 1: Finite size dependence of the pion mass for three values of the quark 
masses. The curves show the exponential fit of eq. (fPfj) with a<i = fixed (see 
text). 
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Figure 2: Finite size dependence of the twist-2 operator matrix element for three 
values of the quark masses. The curves show both the exponential (a 2 = m n fixed) 
and the power law fits (see text). The open symbols denote simulation results with 
lattices of the same spatial extent, but with a shorter time extent (see text). They 
are slightly displaced in a/L for better visibility. 



5 



L/a 


T/a 


N 


12 


36 


3200 


14 


36 


1600 


16 


36 


800 


16 


42 


800 


24 


36 


800 


24 


42 


800 


32 


36 


400 


32 


56 


185 



Table 1: Parameters of our simulation points. N denotes the number of measurements taken 
into account. 





m n (L = oo) 


X 2 /dof 


(x)(L = oo) 


X 2 /dof 


exp: 


re = 0.1340 


0.3533(3) 


1.7 


0.2672(7) 


0.03 




0.1345 


0.2946(3) 


0.6 


0.2579(10) 


0.002 




0.1350 


0.2239(4) 


0.6 


0.2458(24) 


0.12 


exp*: 


re = 0.1340 


0.3533(3) 


2.2 


0.2670(7) 


0.12 




0.1345 


0.2945(3) 


2.0 


0.2576(9) 


0.08 




0.1350 


0.2237(4) 


1.6 


0.2461(15) 


0.09 


pow: 


re = 0.1340 


0.3533(3) 


1.6 


0.2671(6) 


0.03 




0.1345 


0.2946(3) 


0.6 


0.2579(11) 


0.002 




0.1350 


0.2239(4) 


0.6 


0.2453(29) 


0.17 



Table 2: Results for the infinite volume limit of m w and (x) using the exponential fit (exp) of 
eq. (fH|l . fixing <i2 = (exp*) and employing the power law fit of eq. l(T5|l to describe the FSE. 



pow : F(L) = a + ai/L a2 (15) 

The values of m n (L = oo) as obtained from the different fits are summarized in 
table 121 complemented by \ 2 jdoj. In fig. Qwe show the size dependence of m f for 
the three quark masses. The solid curves show the exponential fits keeping a<i = m v 
fixed. Also a power law fit describes the data rather well and leads to a comparable 
number of the inifinite volume pion mass. For the pion mass in the unquenched 
case it is inferred that for very small lattices the FSE are described by a power law 
behaviour jTHj. In this case the source of the FSE is the squeezing of the hadron in 
a small volume. Before the asymptotic (infinite volume) region is reached, there is a 
pre-asymptotic region where the unquenched FSE are described by an exponential 
law P2] indicating that the source of the FSE is the wrapping around the world 
of the virtual pions. Our data obtained in the quenched approximation can not 
differentiate between these two cases. 

In fig.|21we show the finite size dependence of (x) for the three quark masses. As 
expected, the FSE turn out to be stronger for lighter quark masses when keeping 
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Figure 3: Time dependence of the 3-point function for k = 0.1340 coming from the 
following lattice sizes: 32 3 x 56, 32 3 x 36 (squares), 24 3 x 42 and 24 3 x 36 (circles). 



the box size fixed. At present there are no predictions of the form of the FSE for 
a matrix element such as (x). As an ansatz, we have used the same fit functions 
as for the pion mass. Again, the comparison of the matrix element with both fit 
functions does not allow to prefer either of the two fit functions. Thus, carrying 
over the previous arguments for the pion mass to the 3-point function, it is difficult 
to disentangle which is the physical origin of the FSE of the matrix element. 

At this point we want to make a short remark about the importance of the 
physical extent of the plateau in time direction. In fig. El the empty symbols show 
the values of (x) obtained with lattices of the same spatial extent, but with a shorter 
time extent (T" ~ 2.8 fm < T ~ 3.3 — 4.4 fm). Comparing (x) on the lattices 
L/a = 24 and L/a = 32 with smaller time extent, it looks like there would still 
be some FSE. This is, however, not the case and the explanation of this effect has 
a different source. Since it takes around 1.3 fm before all the contributions of the 
excited states go below the accuracy of our computation, a time extent below or 
around 3 fm can be a source of an in-accurate determination of the plateau value, 
given the small number of points available. In fig. |3] we plot the time dependence 
of the 3-point function for k = 0.1340, coming from lattices with different spatial 
extents (L/a = 24, 32) and different time extents (T/a = 36, 42, 56). The two spatial 
extents should be large enough, for all our quark masses, to be free from FSE. Since 
the contributions of the excited states have already died out at around 1.3 fm, even 
if the plateau is short it should give results consistent with the result obtained from 
a longer time extent. From fig. El it is clear that the disagreement between the two 
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Figure 4: Finite size dependence of function of z = m^L. The solid curves 

show the FSE for the pion mass described by formula (|17j) . 

time extents is given by a statistical fluctuation that plays the role of a systematic 
error due to the short time extent. In the longer time extent case too, there is a 
wiggling in the plateau due to statistical fluctuations, but in the short time extent 
case the plateau is not long enough to average away these statistical fluctuations. 
This is also consistent with the fact that for the two different spatial sizes this 
fluctuation goes in opposite directions. A warning we want to give is that if the 
number of points in the plateau is only around 3 — 5 there is the risk to have results 
biased by 2 — 3 a, that obviously are consistent with being a statistical fluctuation. 

It is however quite clear that the FSE for the 3-point functions are much larger 
than for the pion mass. In order to make this statement more quantitative we have 
plotted the relative systematic deviation 

ML) = F{L l7 T F{L 1 °° ) F {L) = m«(L), (x)(L) (16) 
F(L = oo) 

as functions of m^L, where m n is the infinite volume pion mass. Since the FSE for 
the pion mass are very small, it makes no difference to study the dependence of the 
data as functions of m 7Z (L)L or m n (L = oo)L. In fig. 0] and El we see that both the 
data sets scale in a good way as function of m^L. Apart from the L/a = 12 points 
the data sets are in agreement with the formula given in ref . ^§1 120] to leading order 
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Figure 5: Finite size dependence of the twist-2 operator matrix element as a function 
of z = m n L. The solid curves indicate the finite size effects of the pion mass 
described by formula (JT7J) for comparison. 
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of chiral perturbation theory 



8tt 2 \ F n 




) 



2 ATi(m^L) | 2K 1 {y/2m v L) 



3 




2 



4(2tt) 3 / 2 



(rn^L) 3 / 2 (v^ m ^) 3 / 2 



(17) 



where we have interpolated the results for m^/F^ of ref. [18J to our k values. 

A possible explanation for this behaviour is the fact that the validity of the p- 
expansion of chiral perturbation theory, used to determine formula (|17J1 . is limited 
by the constraint F n L 3> 1. For the lattice size L/a = 12 we find instead F n L ~ 1. 

The interesting (and worrisome) result is that while the FSE for the pion mass 
vanish for m n L > 4, the FSE for (x) vanish only for m^L > 5.5. At m n L ~ 4 where 
the FSE for the pion mass are consistent with zero, the 3-point function, due to the 
FSE, is overestimated by around 5%. 

4 Conclusions 

As suggested in [7j, the breaking of the gauge center symmetry Z 3 can be respon- 
sible for larger FSE of the hadron masses in the unquenched case compared to the 
quenched case. This effect, that would lead to power law FSE, is dominant for in- 
termediate volumes. It is also well known that FSE are larger for the nucleon mass 
than for the pion mass [TJ |S] , and moreover, the lattice size where the power law 
behavior disappears in favor of an exponential one is larger for the nucleon than for 
the pion. There is then a danger that the systematic uncertainty coming from FSE 
in the current unquenched determinations of the nucleon matrix elements is severely 
underestimated. At the moment there are no studies of FSE for 3-point functions 
involving twist-2 operators between nucleon states, and to our knowledge, this is 
the first study of FSE of matrix elements involving twist-2 operators. 

We think, it would be important and urgent to have an analytical description of 
the FSE and in particular of the size dependence of the twist-2 matrix elements be- 
tween hadron states, along the line of the computations performed in the framework 
of chiral perturbation theory, in j20j and |21j . 
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